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Introduction 


The Pioneer 10 spacecraft passed within 2.84 Jovian radii (R,) of the planet 

« v jj 

I! 

Jupiter on December 4, 1973. Extensive observations of the Jovian magnetic field and !| 
its interaction with the solar wind plasma were made while the spacecraft was 

I 

within about 100 Rj of the planet. The magnetosphere was found to be severely 
stretched due to the presence of an intense current sheet, which was particularly 
evident during the outbound passage of Pioneer 10 near the dawn terminator jj 

{Smith, et al., 1974). Plots of the angle between the orientation of the outbound 

1 ; 

field and the radius vector from the planet to the spacecraft showed a strong 
tendence for the field to become radial at large distances from the planet (see j! 

Fig. 8, Smith,et al ., 1974). A similar trend has also been seen in both the 

jj 

inbound and outbound Pioneer 11 data (Smith>et al.,1975; Jones, et al.,1975). We jj 

report here some preliminary work cn a mathematical model of the magnetosphere jl 

of Jupiter which is based upon the Pioneer 10 outbound data. A preliminary I; 

model study related to the outbound Pioneer 10 data has also been reported by 
Goertz,et al., (1974). However, we have noted some fundamental conceptual errors | 

in their study and it is also the purpose of this paper to report a correction 
of this earlier analysis. We will also discuss some of the implications of the 
radial field configuration inferred from the Pioneer 10 and 11 data. j 

The Method jj 

Since it is always true that fj 

vH3 =■ 0 ;i 

}j 

one can represent 'B by an expression of the iorm jj 

B = vf X ‘ jj 

where f and g are scalar functions of the coordinates that are sometimes referred 
to as Euler potentials (Euler, 1769; Truesdell, 1954; Stern, 1966). The utility 
of thi c manner of representing B lies in the fact that since B is tangent to the jj 

J' 

intersection of the surfaces f = constant and g = constant, this affords a direct } 


method of evaluating the shape of the lines of force. For example, for a dipole 


in spherical coordinates, we have 

r» M 6 

+ -° * / 



where j> =■ R/g^, and 

3 ' + 

(In cylindrical coordinates,^) will represent the dimensionless component of R 
that is perpendicular to the axis of the dipole.) The f function can be easily 


I 


manipulated into the well-known constant L representation of a dipole field line, 


namely 


s L B 


jAlthough, the law of superposition holds for magnetic fields, this is not generally 


true for the functions f and g, i.e., 


Alternatively, one notes that 

= VX 


so that the vector potential A is related to the f and g functions through 

A = ^ 

For axi -symmetric fields., f is independent of $ and g = . The vector 


potential in spherical coordinates is then 

X = i-CML % 

‘5?V' er * - 2 .VK e ^ 


and in cylindrical coordinates, 

A„ i. = i 


,.~£ 


For axi-symmetric fields, or f functions sharing the same g function, one writes 


A = 

or, alternatively ^ . 

IS 

We write with Bp being the dipole field, 

B =% ^ ^ 2 . 


3 


so that the perturbation field B is given by 

b p » 3 - 2 ? p ■ 

- *&, + 

__ =■ 

B ] represents the axi -symmetric portion of the perturbation field ancft^ the 
component contributing to the spiralling. 


Spherical Polar Coordinate Model 


In spherical coordinates, we have 

T , - -J r - U Jf' 

p si j 




Since there is clearly spiralling of the field (Smith et al . , 1974), we have 
written 

-3 I ^2_ 

- <~p 4- 

where kp represents the spiralling. Because Bj and Bp share the same g function, 

and therefore superposition holds for the respective f functions, we will concentrate 

on these components of the field only. The function F = f^ + f^ will then 

represent meridional plane projections of the measured field. 

In deriving an f-j function, we start with a component of the perturbation 

field whose functional form may be easily deduced from the data. Since the 
radial component of the field decreased and at times reversed, which is con- 
sistent with passage into .3 thin current sheet (Smit^et al., 1974)^ a functional 
form for b. that is consistent with these factors (see also Bird, 1975) is 
l 
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Then 
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q>3= 


we have 


% = 


fa-lib's 

0 ,* fa L- J 


The f function corresponding to the axi -symmetric portion of the perturbation 


field is then given by 

*, - - f ^ 


C-oS & 


l " **#0 y> j 

From the Pioneer 10 outbound data, we find that the constants for f^ are 
approximately 

a = 1.70 
b = 1.10 
A = 7700 
C = 770 

cos e Q = 0:025 

Although fitting the Pioneer 10 outbound data quite well, this F function 

exhibited rather anomalous behavior at high latitudes. Since the expression for 

b 0 can include additional functions of & which are small near e- - T *'/2 (i.e., 

/ \ 

functions of cos fr ) one could write 

t> - _A- + t*os b h (pj A 

f j> a ^*0 ' ' 

Following this lead, alternate functions can be derived. One of several which 


(p *) 


fits the data reasonably well is 
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where the constants are the same as those listed above. The corresponding F 
function is plotted in Figure 1 with M = 4 X 10 5 (Smith, et al., 1974; 1975). 
Although this function exhibits better behavior near the magnetic axis, it 
is still unsatisfactory here. Replacing Ce^“ sin ^ /j>^ by - log cosh(l/cos © Q ) 
produces an which matches the b^ data and is well behaved at 0°, but 
insufficient southward field results because effects due to magnetopause currents 


j 

; f 

, j 
\\ 

|f 

i [ 

J 

i: 

ii 


have not been included in the model. Clearly an additional f function is needed, 
but infinite series techniques will likely be required. 

Neglecting the presence of the magnetopause, the last closed field line 
crosses the magnetic equator at jo c - 440 for the preceding models. Under these 
conditions magnetic field lines originating at higher latitudes would not cross 
the equator and would therefore be considered as being open. 
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The Cylindrical Coordinate Model 


Goertz 4 et al. (1974) have developed a model in cylindrical coordinates. 
For this case, the axi-symmetric portion of the perturbation field is given by 

~ = _ A p 4 1 £ 

*>1 ' 6 ^ 


« 


As before, a functional form for^, that is consistent with the current sheet 
data, etc., Ms a , 

h = * -U.U 

where the current sheet half-width, 0, could be some function of a , i.e. 

= "b „ o/ J 



for a constant angular width sheet, or in general 


Since 


then 



•£.= ~ -^5^ I'h 4 ' c (/0 


The corresponding function for b z is then 

c/ f h 

Note that for z/D^3, log cosh z/D is very nearly z/D - In 2, so that 

-- 


u 

jl 

i 


Plotting^ versus^ for all z/D!>3, and b z versus/) for fixed values of z/D, 
allows one to determine the constants In fj. 6oertz a et al. (1974) have plotted 


b z in this manner and find 
<c*.4 lo\ t =- 2-T7 


so that 


c 0‘) 

+ ' r ) 


However, Goertz et al . (19 74) have plot ted b versus/* , where 

W = ^,/> - u/ 

and 


/b. 


/ 
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As a result, they obtain a power law representation of the component parallel to 
the magnetic equator which lies in the curved surface represented by 
4> 4 kp — CoW^^ayJc 

However, the resulting f will be for such surfaces, does not represent an axi- 
symmetric field, and therefore cannot appropriately be added to the dipole f 
function, which is axi -symmetric. That is, the f functions to be added must 
share the same g function. Goertz, et al . (1974) found that 


and hence 


ra y 


= — V y.icT^ 


- \<> i/Th&ioy j x 


Since they found 

L 




* V> 


over the range 


j j> = 20 to = 80, the corresponding/* dependence for should 


be corrected by the factor z> , or 
I </->■ 72 

K << 

y 

so that 


: 7° 


a = 1.78 


b =-0,01 
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As a check on this, we determined the power law dependence of^ °*y° directly 
and obtained values for a. A, and b of 1.75, 1.0 X 10^, and +0,02 respectively 
for the range ^ = 30 to 80. Combining our results with those of 
Goertz,et al . (1974) we find that f^ is given by 




1.0 

o3s- 

lS 


C & 0 A. -f~ I 


( 30 ^< 80 ) 


where b has been assumed equal to zero, and, based upon one well defined current 
dip* D 0 has been set equal to 1 in units of Jovian radii. The total F function 
representing the axi -symmetric portion of the field is then 


p - 

£ 43 * 






o. 7 r 




A plot of F is shown in Fig. 2. Applying the same conditions as for the 
spherical model, the above model predicts that the last closed field line 
will cross the magnetic equator aty 9 c = 160. Using a * 1 2/3 and A = 7.5 X 10^, 
Goertz, et al . (1974) obtain/). *= 1 50. 


V- 


The Currents 


The current configuration in the magnetosphere can be obtained simply from 
Ampere’s law. Using the field expressions derived from the several fj functions 
one can obtain the configuration of the intense current sheet that exists at the 
magnetic equator as well as the volume currents. The 4 components of the internal 
magnetospheric current system is found in each case to consist of a sheet current 
term plus a volume current term, where the sheet term for the spherical model is 



A st Vi -Q 

AA f 


\€CU 


/<• Prf 
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and for the cylindrical model 



A 

A* 



Although the volume terms are negligible near the magnetic equator, they dominate 

near the magnetic polar axis. This is clearly an artifact of each model which 

k 

will disappear when proper account of the magnetopause currents and ''the inner 

, , it Asc 

cutoff radius of the current srieet Aare included. 

Mathematically terminating the field at the magnetopause allows one to 
solve for the magnetopause currents, and the corresponding boundary field direction 



if 
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can be determined for the several models and compared with the data. For 
example, a function, ttj 1 ) , which can terminate the field arbitrarily abruptly 

1 + Q~ %o) 


is 


W " 


so that the terminated field, B', is given by 

s' ■ 

Here J> o is the radial distance to the magnetopause (as a first approximation 

we assume the- magnetopause boundary to be spherical) and d relates to the thickness 

of the boundary. In principle, B should be the total field. However, we still 

neglect the component of the perturbation field. 

The above function terminates the preceding azimuthal currents at J> ~j> o 

and addition provides the magnetopause currents,- i .e. , for the spherical 
model 




*4- 






X ('~ fco) 


>0 


and for the cylindrical model (here J> - l/ + y2 , normalized) 
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_ ot 






' f° K 




we find that B^ for the spherical model is positive at all values of 0 so that 
the corresponding magnetopause current is clockwise, as viewed from the magnetic 

pole, at all latitudes. Hence, just prior to the magnetopause boundary the pred- 
icted field direction is southward,- as is observed by both Pioneers 10 and 11 
(Smith, eL al . , 1 974 ;1 975 ) . On the other hand, the bracketed term contained in 


11 


the magnetopause expression for the cylindrical model becomes negative at 

magnetic latitudes greater than about 20° so that the direction of the magnetopause 

current flow reverses from a clockwise direction at lower latitudes to a 

counter clockwise direction at higher latitudes. The corresponding field 

just inside the boundary is predicted to point northward at latitudes greater 
_ h 

than 20° and sou^ard at lower latitudes. Such a prediction appears to be 
in disagreement with the data although it is interesting to note that prior to 
the outbound magnetopause crossings by Pioneer 11 there were intervals 
approaching 8 hours during which the magnetospheric fields pointed north of 
radial by roughly 20°. However, these are likely transient features related 
to a (possibly) northward component of the solar wind flow velocity. 


Discussion 


Since the functions plotted in Figs. 1 and 2 were derived from the 
Pioneer 10 outbound data, they were developed from data taken within about 
20° of the magnetic equator and over the radial range 20 <£ 80, and 

qualitatively represent the magnotospheric field configuration in meridional 
planes that lie near the dawn terminator. However, one notes that the 
models also qualitatively fit thePPioneer 11 data quites well which extends 
the latitude range of the functions to perhaps 40° (the Pioneer 11 inbound data 
is qualitatively very similar to the Pioneer 10 outbound data) and to about 
40° of the dawn meridian. As is evident from the figures, both models 

should be considered unreliable at latitudes greater than about 45°. 

A basic difference between the two models is the fact that one is for a 
constant angular width current sheet (the spherical coordinate model, Fig. 1) 
while the other is for a constant thickness current sheet (the cylindrical 
coordinate model. Fig. 2). Likely the actual case lies somewhere between these 
two current sheet configurations. Both permit the current sheet to exist to 

• v > 

the center of the planet although. it ‘ must be cut off at some 

inner radiuy>,> 2, since one would not expect the sheet to exist within the 
centrifugal -gravitational “balance distance of several radii. 

Another basic difference involves the direction of flow of the magnetopause 
currents and the corresponding direction of the magnetopause field. Predictions 
based upon the spherical model are more consistent with the measurements. 

There are also a number of factors regarding the constants derived for 
the models that should be mentioned. For example, in the case of the cylindrical 
coordinate model, the value of c in fj is determined from b z versus at constant 
z/D, but this requires a knowledge of D. The evaluation of the constants a. A, 
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and b depends c 'tically upois the accurate determination of the actual power law 

dependence of D on . The current sheet half-width is one of the most uncertain 

parameters and its dependence upon y> is particularly difficult to determine 
directly from the plots of B versus ^ . The constants contained in the 
expression for fj are self-consistent and the model appears to establish the 
yd independence of D. Similar comments can be made regarding the constant 
angular width model as well. Unfortunately, a brief study of the variation of the 

widths of the field dips has not shed much light on this crucial point except 

that the data tend to favor the cor-tant angular width model. 

In the study by Goertz et al . (1974), the determination of a and D (their 
b 0 )from a plot of b s instead of b^ causes the resulting function that is to 
represent the shape of the field in meridional planes to be a mixture of f 
functions requiring different g functions. On the other hand, in our determination 
of these constants for the cylindrical model, we have used only f functions that 
have the same g function. Because the spiralling of the field was not excessive, 
the disagreement with the results of Goertz et al. (1974) is not great, and a 
comparison of the plots of the field lines shows them to be quite similar. 

Any interpretation regarding the value ofy> c (where By= 0) that is derived 
from the models should be viewed with caution since the model fits the data only 
out to about 80 or 85 Rj, and hence these cutoff radii should be considered as 
possible artifacts of the models. An artifact of this kind is meaningless 
because the magnetopause currents have been neglected in the derivation of the 
f functions. As noted earlier, it is tempting to assume that field lines leaving 
the planet at higher magnetic latitudes than those related ti o yj c are open field 
lines and that they merge with the interplanetary field. But the data show that 
the field lines are southward at the magnetopause, suggesting that they are closed 
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by the magnetopause currents. In the sense of field lines and particle trapping, 
these lines clearly will not have trapped particles on them. The last closed 
field line which could contain trapped particles should be the one which crosses 
the equator just prior to the magnetopause boundary. 

The particles in the intense equatorial current sheet likely result from 
plasma flow due to the combined action of a Jovian "polar wind," much like 
that postulated for Earth (Banks and Holzer, 1969), plus the strong centrifugal 
force caused by the large size and rapid rotation of the magnetosphere. 

The balance of pressures at the magnetopause likely must include that exerted 
by a radial flow of polar wind ions moving parallel to the essentially radial 
field lines in the magnetosphere. Perhaps such t plcsnw flow also provides 
a significant stabilizing influence for the large seals magnetosphere configuration 
reported here, since one would otherwise expect the seW u.d to blow the high 
latitude field lines back into the tail because of •or* relatively weak magnetic 
pressure exerted at the magnetopause (Smith, et a:., ; '?4). A study of the 
Pioneer 11 outbound data will provide some important information in this regard, 

although much higher latitude data are clearly needed. 

Further studies of the magnetosphere will likely require the use of 
perturbation techniques (Stern, 1967) in order to obtain more well behaved 
functions at high latitudes and to allow for a non-spherical magnetopause 
boundary. Other studies being conducted at the present time will merge models 
developqjfor the range \ 6 with the magnetospheric models reported here. 

In tliis regard, magnetospheric studies establish reasonable estimates of the 
magneLospheric current systems and detailed attempts at merging the two 
programs will establish, among other things, the inner cutoff radius of the 
current sheet. 

This research was supported in part by the National Aeronautics and Space 
Administration under NASA-Ames contract NAS2-7358. 
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Figure Captions 


t 

Meridional plane representation of the magnetospheric field as developed 

i 

in a spherical coordinate representation. The solid and dashed arrows 

represent the average field direction measured at several points along 

the outbound trajectories of Pioneer 10 and 11, respectively. The 

shaded half-angular width portion near the magnetic equator represents 

a portion of the equatorial current sheet configuration assumed in this 

model, and a first order approximation to the magnetopause boundary is 

also indicated. Regions over which the function is reliable are % 
indicated in the text. The curves leave the j> - 2 sphere at equally 

spaced angular intervals. 

Meridional plane representation of the magnetospheric field as developed 
in a cylindrical coordinate representation. The solid and dashed arrows 
represent the average field direction measured at several points along 
the outbound trajectories of Pioneer 10 and 11 respectively. The 
shaded half-width portion near the magnetic equator represents a portion 
of the equatorial current sheet configuration assumed in this model, and 
a first order approximation to the magnetopause boundary is also indicated. 
Regions over which the function is reliable are indicated in the text. 

The curves leave the = 2 sphere at equally spaced angular intervals. 
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